How To Generate As Many Somos-Like Miracles as You Wish 

Shalosh B. EKHAD and Doron ZEILBERGER^ 

Somos Sequences 

In the late 1980s, Michael Somos came up with an amazing set of conjectures, featuring non-linear 
recurrences that always yield integers. The most famous one was Somos-4: [S]http : //oeis . org/A006720. 
^ . Define a sequence {a(n)}^]^ by a(l) = 1 , a(2) = 1 , o(3) = 1 , a(4) = 1, and for n > 5, 

O' a{n - l)a{n - "i) + ain - 2f 

(N ■ a(n) = . 

a[n — 4) 

H ■ 

^ ' Prove that a{n) is an integer for every n > 1. This is indeed amazing, since, starting at n = 9, you 

divide by an integer larger than 1, and a-priori you are only guaranteed to get rational numbers, 
but surprise-surprise, you always (seem) to get integers. 



(N 

o 
u 



This took quite some effort to prove rigorously. Of course, confirming it for n < 300, say, already 
constitutes a (very convincing!) empirical proof, but mathematicians do not care about truth per- 
se, they only care about playing their (sometimes fun but often very dull) artificial game called 
proving ( "rigorously" , using logical deduction) . 

And indeed, such a mathematical (elegant!) proof was given by Janice Malouf[M]. Other people 
gave other proofs, see the references given at [S] http://oeis.org/A006720. 

Of course, Michael Somos is a honest person. He first encountered this sequence in some context. 



I and then made his conjecture. But suppose that you are dishonest, and want to amaze your friends 

' (or challenge your foes) with such Somos Phenomena, coming up with non-linear recurrences where 



you know beforehand, that you are guaranteed to get integers, and even have a fully rigorous proof 
' of that fact. 



Then you must follow the advise of Carl Gustav Jacob Jacobi who taught us that in order to make 
life simpler. 



, ^^man muss immer umkehren" , 

in other words, work backwards! Instead of trying to play darts, by placing the target on the wall, 
and then aiming from a distance of twenty meters, and most probably missing it, one first throws 
the dart, and then draws the bull's eye! 

Consider the linear recurrence sequence defined by 

/(l) = l,/(2) = 3,/(n) = 3/(n-l)-/(n-2) (n > 3) . 
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NesSomos downloadable from Zeilberger's website. Supported in part by the NSF. 
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Obviously (by induction on n) they are all integers. 



Consider, on the other hand, the foUowing sequence 

6(l) = l,6(2) = 3,6(n) = ^^^i^-ll^ (n > 3) , 
now it is not at all immediately obvious that the 6(n) are always integers. 

Suppose that one asked the question "prove that the b{n) are always integers" . The best way to 
prove it would be to generate the first few terms, go to Neil Sloane's OEIS[S], and find out that 
6(n) = A001519(ra). In other words it is the same as the above sequence f{n), at least for the first 
30 terms. 

So far it is only a "conjecture", but if you know about the C-finite ansatz ([K][KP][Z1][Z2]), then 
knowing that f{n) = h{n) for 1 < n < 4 implies that it is true for all n > 1, thereby yielding an 
'empirical', yet fully rigorous, proof! 

It is easy to see that C-finite sequences, i.e. sequences that satisfy linear recurrence equations with 
constant coefficients, also satisfy (many!) non-linear recurrences. But more surprisingly, the same 
is true for any subsequence where the places are either polynomials (e.g. n^) or exponential 
(e.g. 2"'), or even exponential polynomial (e.g. n^2" + n'*3"). For example, if F{n) are the 
Fibonacci numbers defined by F{0) = 0,F{1) = l,F{n) = F{n - 1) + F{n - 2)(n > 2), then the 
sequence 

a(n) := ^(71^2" + n=^5" + 1) 

satisfies some (very complicated!) non-linear recurrence equation, i.e. there exists a positive integer 
r (the order) and a polynomial P{xo, . .. ,Xr) such that 

P{a{n),...,a{n + r)) = {for all n > 1). 

Furthermore, by possibly increasing r, one can demand that the degree in Xr is 1, so one can express 
o(n + r) as a rational function of a(n), . . . , a(n + r — 1) and get a genuine non-linear recurrence 
that very surprisingly (if you don't know how it was found!) yields integers. 

Of course the same is true for any sequence that satisfies a linear recurrence equation with constant 
coefficients. 

So let's state the main theorem in its intimidating full generality. 

Main Theorem: Let /(n) be a sequence satisfying a linear recurrence equation with constant 
coefficients 

d 

f{n) = ^Cif{n - i) , 
1=1 

for some constants ci, . . . , cj, and with some given initial conditions /(I), • • • , f{d). Let p{n) be of 
the form, for some specific integer m, 

p(n) := POL(n,2",3",...,m") , 
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where POL is some polynomial of m variables. Then there exists an integer r and a polynomial 

P(a:o, . . . , Xr), of degree 1 in Xr such that the subsequence 
a{n) := fip{n)) satisfies the non-linear recurrence 

P{a{n),...,a{n + r)) = (n > 1) . 



The formal proof of this theorem would be very boring, but its informal proof, at least generically, 
is not too bad. Using Binet's formula, one can express f{n) "explicitly" as a linear combination of 
a" (and possibly na'^,n^a^ etc. in case of multiple roots) where the a's are roots of the equation 

d 
i=l 

By using the initial conditions, one can find the coefficients in that linear combination of the a^'s, 
that are also certain expressions in the a's. Then, since q;("'+^)^ = a'^^ {a" )'^a, a^"^^ = (a^" )^ etc. 
one can introduce auxiliary variables, {a", a" ,0" , . . . ,■■■}, (finitely many!), and recall 

that the a's all satisfy the algebraic equation above (and enjoy nice relations amongst themselves, 
e.g. that their sum is ci and their product is {—iy~^Cd)- 

Now one can express f{p{n)) in terms of these auxiliary variables, and by elimination (in principle 
using the Buchberger algorithm), one can get, for some r, a pure polynomial relation, of the desired 
form P(a(n), . . . , a(n + r)) = 0, and even demand that the degree in a{n + r) is one, so that one 
can express a{n + r) as a rational function of a(n), . . . , a{n + r — 1). 

In practice, it is better to use undetermined coefficients. Fixing r and d, write a generic polynomial 
P{xo, . . . , Xr) of total degree d (and degree 1 in Xr), plug-in the expressions in terms of the auxiliary 
variables, get a huge polynomial in these variables, set up all the coefficients (w.r.t. the auxiliary 
variables) of the numerator to 0, and solve the large system of linear equations where the unknowns 
are the "undetermined coefficients" , thereby hopefully making them determined. If the only solution 
is the trivial zero solution, don't give up! Just raise d and/or r. 

Now if we take the coefficients of the recurrence ci, . . . , to be integers, and the initial conditions 
/(I), . . . , f{d) to be integers as well, then the sequence /(n) obviously consists of only integers. If, 
in addition, the polynomial POL has all positive integer coefficients, then we are guaranteed, a 
priori, that the terms of the subsequence a(n) := f{p{n)) are all integers! 

For example, if our sequence is /(n) = i-e. /(I) = 1 , /(2) = 3 , 

/(n) = 3/(n — 1) — /(ra — 2) (n > 3), considered above, then the Binet formula is 

„, , a" -a-" 
/(«) 



a — a ^ 



where a is a root of the quadratic equation a — 3a -|- 1 = 0. Taking p{n) = n we get 

a(n) 



a — a ^ 



3 



Now 

a{n + 1) = 



a — a ^ a — a ^ 



an + 2 = 3^ = ^ . 

a — a a — a 

So you have expressions for a(n),a{n + 1), a{n + 2), as rational functions in the auxiliary variables 
A := a", B := a" (and a, but a satisfies — 3a + 1 = 0). So generically we should have a 
second-order polynomial relation P{a{n),a{n + l),a{n + 2)) = 0. But it is not of degree one in 
a{n + 2). If we allow higher-order recurrences, it turns out that we have a fairly nice non-linear 
recurrence of order 5. Here is one such example. 

Prop. Let a(n) be defined by 

a(l) = 1 , a(2) = 21 , a(3) = 2584 , a(4) = 2178309 , a(5) = 12586269025 

and for n > 6, by the recurrence 

a(n) = ( 2303 a{n - 4)a(n - 3)a(n - 1) + 2255 a{n - Sfa{n - 2) + 329a(n - 4)a(ra - 1)^ 

- 15792 a{n - 4)a(n - 2f + 329 a{n - 4)a(n - 3)^ - 2303 a{n - Afa{n - 2) + 441 a{n - 2) 
— a(n— 2)^^-2961 a(n— 4) — a(n— 5)a(n— 2)a(n— 1) + 329a(n— 5)a(n— 3)o(n— 2) )/(48a(n— 4)o(n— 2) ) , 
then a(n) are all integers! 

For many more such propositions, see: 

http : //www .math . rutgers . edu/~zeilberg/tokhniot/oNesSomos3 . 

The much sparser subsequence, a(n) := /(2") satisfies a much simpler (second-order!) non-linear 
recurrence: 

a(n- 1) f4-F5a(n- 1)^) 

a(n) = ^ 5 ^ 

2 + 5a(n-2)^ 

subject to the initial conditions 

a(l)=3 , a(2) = 21 . 

The Maple package NesSomos 

Everything is implemented in the Maple package NesSomos (written by DZ) available directly from 
http : //www . math . rutgers . edu/~zeilberg/tokhniot/NesSomos , 

and that is linked to from the "front" of this article 

http://www.math.rutgers.edu/~zeilberg/mamarim/mamarimhtml/somos.html , where one can 
also find three long webbbooks generated by SBE. 
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To get a list of the main procedures, type ezraO ;, and for help with a specific procedure, type 
ezra(ProcedureNaiiie) ;. Let's just list the more important procedures. 

• FindREL(LPols,Vars,x,d) is a very general procedure that inputs a list of rational functions 
LPols in the list of variables Vars, a symbol, x, and a positive integer d, and outputs a polynomial 
of degree < d, let's call it P(a;[l], . . . ,x[r\) such that 

P{LPols[l],...,LPols[r]) = Q . 

For example, 

FindREL{\p^ - q^,2pq,p^ +q\\p,q],x,2)- 

yields 

2 2 2 n 

If nothing is found, then it returns FAIL and one has to increase d. Of course in order to guarantee 
that there is such a good d, we must insist that the length of LPols is larger than the length of 
Vars. 

• FindlnREL (LPols, Vars, x,d) is exactly as above, but, in addition, the output should have degree 
one in x[r\ (where r is the length of LPols)). This is needed for the application to generating Somos- 
like sequences. 

• FindSMpG(c,p,n,Maxr,d,a) ; inputs a positive integer c, a polynomial expression p in the discrete 
variable n, positive integers Maxr and d, and a symbol a, and outputs a non-linear recurrence, of 
order Maxr and degree < d satisfied by the sequence a{n) := f{p{n)), where f{n) is the sequence 
satisfying /(n) = c/(n — 1) — /(n — 2) and /(O) = 0, /(I) = 1. If it fails it returns FAIL, and the 
human has to try again with a larger d and/or Marx. The first amazing example above was gotten 
from 

FindSMpG (3 , n**2 , n , 5 , 6 , a) ; 

• FindSMeG(c,e,n,Maxr,d,a) ; inputs a positive integer c, a positive integer e, positive integers 
Maxr and d and a symbol a, and outputs a non-linear recurrence of order Maxr and degree < d 
satisfied by the sequence fie"-), where /(n) is the sequence satisfying /(n) = c/(n — 1) — f{n — 2) 
and /(O) = 0, /(I) = 1. If it fails it returns FAIL. 

The second amazing example above was gotten from: 

FindSMeG(3,2,n,3,3,a) ; 

Webbbooks 

For many amazing Somos-like miracles see the three webbooks 
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• http : //www .math . rutgers . edu/~zeilberg/tokhiiiot/ oNesSomosl 

• http : //www .math . rutgers . edu/~zeilberg/tokhniot/oNesSomos3 

• http : //www .math . rutgers . edu/~zeilberg/tokhiiiot/ oNesSomos4 

Enjoy! 

Using NesSomos you can find many more Somos-like miracles on your own, impressing your friends 
and challenging your enemies. 

Laurent Phenomenon 

Analogous things can be done to generate sequences that have the so-called Laurent phenomenon, 
and even, a priori, the much stronger property of positivity. 

Encore 

By looking at the output of 

http://www.math.rutgers.edu/~zeilberg/tokhniot/oNesSomosl , we see a pattern that can 
be summarized as follows. 

Proposition: Let c be a positive integer (a formal symbol), and define a sequence by 

a(0) = 1 , a(l) = c , 

and for n > 2, 

a(n - 1) (4 + (c - 2)(c + 2)a(n - 1)^ ) 
2 + (c-2)(c + 2)a(n-2)2 

Then for all n > 1, a(n) are integers (polynomials in c with integer coefficients). 
Proof: Consider the sequence defined by /(O) = 0, /(I) = 1 and for n > 2 by 

/(n) = c/(n - 1) - /(n - 2) . 

Define 6(n) := /(2"). It is routine to prove that 

b{n - 1)(4 + (c - 2)(c + 2)6(n - 1)^) 



6(n) = 



2 + (c-2)(c + 2)fe(n-2)2 



and, of course, 6(0) = 1 and 5(1) = c. Hence, by induction, b{n) = a{n) for all n > 1. Since the 
/(n) are obviously integers (polynomials in c with integer coefficients), so are a(n). □ 
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